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ABSTRACT 

The X and 7-ray observations of the source GRB 031203 by INTEGRAL are interpreted within our theo- 
retical model. In addition to a complete space-time parametrization of the GRB, we specifically assume that 
the afterglow emission originates from a thermal spectrum in the co-moving frame of the expanding baryonic 
matter shell. By determining the two free parameters of the model and estimating the density and filamentary 
structure of the ISM, we reproduce the observed luminosity in the 20-200 keV energy band. As in previous 
sources, the prompt radiation is shown to coincide with the peak of the afterglow and the luminosity substruc- 
ture are shown to originate in the filamentary structure of the ISM. We predict a clear hard-to-soft behavior 
in the instantaneous spectra. The time-integrated spectrum over 20 seconds observed by INTEGRAL is well 
fitted. Despite this source has been considered "unusual", it appears to us a normal low energetic GRB. 
Subject headings: gamma rays: bursts — gamma rays: observations — radiation mechanisms: thermal 



1. INTRODUCTION 

GRB 031203 wa s observed by IBIS on b oard of the IN- 
TEGR AL satellite llMereghetti & G6tsll2003l). as well as by 
XMM iTWatson et alJ 120041) and Chandra ISoderberg et all 
2004) in the 2- 10 keV band, and by VLT dSoderberg et alJ 
2004) in the radio band. It appears as a typical long 
burst (Sazonov. Lutovinov & Sunvaev 2004a), with a sim- 
ple profile and a duration of w 40 s. The burst fluence 
in the 20-200 keV band is (2.0 ± 0.4) x 10" 6 erg/cm 2 
ISazonov. Lutovinov & Sun vaevl l20 04al). an d the measured 
redshift is z = 0.106 (Prochaskaet al. 2004). We analyze in 
the following the gamma-ray signal received by INTEGRAL. 
The observations in o ther wavelengths, in analogy w ith the 
case of GRB 980425 JPian et alJl2000t IRuffini etafl l2004b). 
could be related to the supernova event, as also suggested by 
ISoderberg et alJd2004l) . and they will be examined elsewhere. 

The INTEGRAL obser vations find a direct explanation in 
our theoretical model (see Ruffini et al. 2001a b, 2003, 2005c; 
iBianco & R uffini 2005a b| and references therein). We de- 
termine the values of the two free parameters which char- 
acterize our model: the total energy stored in the Dyado- 
sphere Ed ya and the mass of the baryons left by the collapse 
MflC 2 = BEdya- We follow the expansion of the pulse, com- 
posed by the electron-positron plasma initially created by the 
vacuum polarization process in the Dyadosphere. The plasma 
self-propels itself outward and engulfs the baryonic remnant 
left over by the collapse of the progenitor star. As such pulse 
reaches transpa rency, the Proper Gamma-Ra y Burst (P-GRB) 
is emitted (IRuffini et alJll999l 120001 l2001bl) . The remaining 
accelerated baryons, interacting with the interstellar medium 
(ISM), produce the afterglow emission. The ISM is described 
by the two additional parameters of the theory: the average 
particle number density < uism > and the ratio < 1Z > be- 
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tween t he effective emitting area and the total area of the pulse 
( IRuffini et alJl2004al). which take into account the ISM fila- 
mentary structure iRuffini et all20 05a). 

We reproduce correctly in several GRBs and in this spe- 
cific case (see e.g. Fig. Q t he observed time variability of 
the prompt emission (see e.g. Ruffin i et alJ2002tl2003l l2005c. 
and references therein). The radiation produced by the inter- 
action of the accelerated baryons with the ISM agrees with 
observations both for intensity and time structure. 

The progress in reproducing the X and 7-ray emission as 
originating from a thermal spect rum in the co-moving frame 
of the burst (Ruffini et al. 2004a) leads to the characterization 
of the instantaneous spectral properties which are shown to 
drift from hard to soft during the evolution of the system. The 
convolution of these instantaneous spectra over the observa- 
tional time scale is in very good agreement with the observed 
power-law spectral shape. 

As shown in previous cases (see Ruffi ni et alJ2003l l2005b). 
also for GRB 031203, using the correct equations of motion, 
there is no need to introduce a collimated emission to fit the 
afterglow observations (see also Soderberg et al. 2004J who 
find this same conclusion starting from different considera- 
tions). 

2. THE INITIAL CONDITIONS 

The best fit of the observational data leads to a total energy 
of the Dyadosphere Ed ya = 1 .85 x 10 50 erg. Assuming a black 
hole mass M = 10M Q , we then have a black hole charge to 
mass ratio £ = 6.8 x 10~ 3 ; the plasma is created between the 
radii r\ = 2.95 x 10 6 cm and = 2.81 x 10 7 cm with an initial 
temperature of 1.52 MeV and a total number of pairs N e + e - = 
2.98 x 10 55 . The amount of baryonic matter in the remnant is 
B = 7.4xl0" 3 . 

After the transparency point and the P-GRB emission, the 
initial Lorentz gamma factor of the accelerated baryons is 
7 = 132.8 at an arrival time at the detector t d a = 8.14 x 10" 3 
s and a distance from the Black Hole r = 6.02 x 10 12 cm. 
This corresponds to an apparent superluminal velocity along 
the line of sight of 2.5 x 10 4 c. The ISM parameters are: 
< n ism >= 0.3 particle/cm 3 and <TZ>= 7.81 x 10" 9 . 

3. THE GRB LUMINOSITY IN FIXED ENERGY BANDS 
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FIG. 1. — Theoretically simulated light curve of the GRB 031203 prompt 
emission in the 20-200 keV energy band (solid red line) is compared with the 
observed data (green points) from Sazonov, Lutovinov & Sunvaev 1 2004a). 
The vertical bold red line indicates the time position of P-GRB. 




FIG. 2. — Five different theoretically predicted instantaneous photon num- 
ber spectrum N(E) for t/ t =2, 6, 10, 14, 18 s are here represented (colored 
curves) together with their own temporal convolution (black bold curve). The 
shapes of the instantaneous spectra are not blackbodies due to the spatial con- 
volution over the EQTS (see text). 



The aim of our model is to derive from first principles 
both the luminosity in selected energy bands and the time re- 
solved/integrated spectra. The luminosity in selected energy 
bands is evaluated integr ating over the equit emporal surfaces 
(EQTSs, see lBianco & Ruffinill2004l200 5a) the energy den- 
sity released in the interaction of the accelerated baryons with 
the ISM measured in the co-moving frame, duly boosted in 
the observer frame. The radiation viewed in the co-moving 
frame of the accelerated baryonic matter is assumed to have a 
thermal spectrum and to be produced by the interaction of the 
ISM with the front of the expanding baryonic shell. 

In order to evaluate the contributions in the band [v\ , v{\ we 
have to multiply the bolometric luminosity with an "effective 
weight" W(v\,V2 1 T arr ), where T arr is the observed tempera- 
ture. W{v\,V2 1 T arr ) is given by the ratio of the integral over 
the given energy band of a Planckian distributio n at tem pera- 
ture T arr to the total integral aT* rr jRuffini et alJl2004al) . The 
resulting expression for the emitted luminosity is 



dt^dil 



EQTS 



^ v cost? A" 4 ^W{v u v 2 J arr )dY<, 



(1) 

where Ae = AEj„ t /V is the energy density released in the in- 
teraction of the accelerated baryons with the ISM measured 
in the co-moving frame, A = 7(1 -(v/c)cosi3) is the Doppler 
factor, dY, is the surface element at detector a rrival time t d n 
on whi ch the integration is performed (details in Ruffi ni et al.l 
120041 

4. THE GRB 031203 "PROMPT EMISSION" 

In order to compare our theoretical prediction with the ob- 
servations, it is important to notice that there is a shift between 
the initial time of the GRB event and the moment in which the 
satellite instrument has been triggered. In fact, in our model 
the GRB emission starts at the transparency point when the 
P-GRB is emitted. If the P-GRB is under the threshold of the 
instrument, the trigger starts a few seconds later with respect 
to the real beginning of the event. Therefore it is crucial, in the 
theoretical analysis, to estimate and take into due account this 
time delay. In the present case it results in At d = 3.5 s (see the 
bold solid line in Fig. 0. In what follows, the detector arrival 
time is referred to the onset of the instrument. 



The structure of the prompt emission of GRB 031203, 
which is a single peak with a slow decay, is reproduced as- 
suming an ISM which has not a constant density but presents 
several density spikes with < hism >= 0. 16 particle/cm 3 . Such 
density spikes corresponding to the main peak are modeled 
as three spherical shells with width A and density contrast 
An/n: we adopted for the first peak A = 3.0 x 10 15 cm 
and An/n = 8, for the second peak A = 1.0 x 10 15 cm and 
An/n =1.5 and for the third one A = 7.0 x 10 14 cm and 
An/n =1. To describe the details of the ISM filamentary 
structure we would require an intensity vs. time information 
with an arbitrarily high resolving power. With the finite res- 
olution of the INTEGRAL instrument, we can only describe 
the average density distribution compatible with the given ac- 
curacy. Only structures at scales of 10 15 cm can be identi- 
fied. Smaller structures would need a stronger signal and/or a 
smaller time resolution of the detector. The three clouds here 
considered are necessary and sufficient to reproduce the ob- 
served light curve: a smaller number would not fit the data, 
while a larger number is unnecessary and would be indeter- 
minable. 

The result (see Fig. shows a good a greement with 
the ligh t curve reported by Sazo nov. Lutovinov & Su nvaev 
(2004ai), and it provides a further evidence for the possibil- 
ity of reproducing light curves with a compl ex time variabil- 
ity thr ough ISM inhomogeneities jRuffini et al J 120021 120031 
2005c) (s ee also the analysis of the prompt emission of GRB 
991216 in lRuffini et alJ2002l) . 

5. THE GRB 031203 INSTANTANEOUS SPECTRUM 

As outlined in section [5] in addition to the the luminosity 
in fixed energy bands we can derive also the instantaneous 
photon number spectrum N(E). In Fig. |2]are shown samples 
of time-resolved spectra for five different values of the arrival 
time which cover the whole duration of the event. 

It is manifest from this picture that, although the spectrum 
in the co-moving frame of the expanding pulse is thermal, the 
shape of the final spectrum in the labor atory frame is clearly 
non thermal. In fact, as explained in Ruffini et al. (2004a), 
each single instantaneous spectrum is the result of an integra- 
tion of hundreds of thermal spectra over the corresponding 
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FIG. 3. — The energy of the peak of the instantaneous photon number 
spectrum N(E) is here represented as a function of the arrival time during the 
"prompt emission" phase. The clear hard to soft behavior is shown. 



EQTS. This calculation produces a non thermal instantaneous 
spectrum in the observer frame (see Fig. |3. 

Another distinguishing feature of the GRBs spectra which 
is also present in these instantaneous spectra, as shown in 
Fig. 13 is the hard to soft transition during the evolution of 
the event JCrider et alJl997tlPirarJl999HFrontera et all2000t 
iGhirlanda. Celotti & Ghiselhnil2002l) . In fact the peak of the 
energy distributions E p drift monotonically to softer frequen- 
cies with time (see Fig.|3jl. This feature explains t he change in 
the power-law low energy spectral index a ( Ban d et all!993l) 
which at the beginning of the prompt emission of the burst 
(t% = 2 s) is a = 0.75, and progressively decreases for later 
times (see Fig. 0. In this way the link between E p and a iden- 
tified by Crider et al. ( 1997) is explicitly shown. This theoreti- 
cally predicted evolution of the spectral index during the event 
unfortunately cannot be detected in this particular burst by 
INTEGRAL because of the not sufficient quality of the data 
(poor photon statistics, see Sazo nov. Lutovinov & Sunvaevl 
120041). 

6. THE GRB 031203 TIME-INTEGRATED SPECTRUM AND THE 
COMPARISON WITH THE OBSERVED DATA 

The time-integrated observed GRB spectra show a 
clear power-law behavior. Within a different frame- 
work Shakura, Sunyaev and Zel'dovich (see e.g. 
Pozdniakov, Sobol & Sunvaev 1983, and references therein) 
argued that it is possible to obtain such power-law spectra 
from a convolution of many non power-law instantaneous 
spectra evolving i n time. This result was recalled and ap- 
plied to GRBs by Blinnikov. Kozvreva & Panchenko ( 1999) 
assuming for the instantaneous spectra a thermal shape with 
a temperature changing with time. They showed that the 
integration of such energy distributions over the observation 
time gives a typical power-law shape possibly consistent with 
GRB spectra. 

Our specific quantitativ e model is more complicated than 
the one considered by Blinnikov. Kozvreva & Panchenko 
( 1999): as pointed out in section [5j the instantaneous spec- 
trum here is not a black body. Each instantaneous spec- 
trum i s obtained by an integration o ver the corresponding 
EQTS (Bia nco & Ruffinill2004t 12005 al) : it is itself a convolu- 
tion, weighted by appropriate Lorentz and Doppler factors, of 
~ 10 6 thermal spectra with variable temperature. Therefore, 
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FIG. 4. — Three theoretically predicted time-integrated photon number 
spectra N(E) are here represented for < t£ < 5 s, 5 < < 10 s and 
10 < < 20 s (colored curves). The hard to soft behavior presented in Fig. 
l3*lis confirmed. Moreover, the theoretically predicted time-integrated photon 
number spectrum N(E) corresponding to the first 20 s of the "prompt emis- 
sion" (black bold c urve) is comp ared with the data observed by INTEGRAL 
(green points, see Sazonov, Lutovinov & Sunvaev 2004a b). This curve is 
obtained as a convolution of 108 instantaneous spectra, which are enough to 
get a good agreement with the observed data. 



the time-integrated spectra are not plain convolutions of ther- 
mal spectra: they are convolutions of convolutions of thermal 
spectra (see Fig. |2ji. 

The simple power-law shape of the integrated spectrum is 
more evident if we sum tens of instantaneous spectra, as in 
Fig. @] In this case we divided the prompt emission in three 
different time interval, and for each one we integrated on time 
the energy distribution. The resulting three time-integrated 
spectra have a clear non-thermal behavior, and still present 
the characteristic hard to soft transition. 

Finally, we integrated the photon number spectrum N(E) 
over the whole duration of the prompt event (see again Fig. 
0}: in this way we obtain a typical non-thermal power-law 
spectrum which res ults to be in good agreement w ith the IN- 
TEGRAL data (see lSazonov. Lutovinov & Sunvae\l2004albl) 
and gives a clear evidence of the possibility that the observed 
GRBs spectra are originated from a thermal emission. 

7. CONCLUSIONS 

We show how, applying our model to the GRB 031203, we 
are able to predict the whole dynamic of the process which 
originates the GRB emission fixing univocally the two free 
parameters of the model, Ej ya and B. Moreover, it is possible 
to obtain the exact temporal structure of the prompt emission 
taking into account the effective ISM filamentary structure. 

The important point we like to emphasize is that we can 
get both the luminosity emitted in a fixed energy band and 
the photon number spectrum starting from the hypothesis 
that the radiation emitted in the GRB process is thermal in 
the co-moving frame of the expanding pulse. It has been 
clearly shown that, after the correct space-time transforma- 
tions, both the time-resolved and the time-integrated spectra 
in the observer frame strongly differ from a Planckian distri- 
bution and have a power-law shape, although they originate 
from strongly time-varying thermal spectra in the co-moving 
frame. We obtain a good agreement of our prediction with 
the photon number spectrum observed by INTEGRAL and, 
in addition, we predict a specific hard-to-soft behavior in the 
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instantaneous spectra. Due to the possibility of reaching a pre- 
cise identification of the emission process in GRB afterglows 
by the observations of the instantaneous spectra, it is hoped 
that further missions with larger collecting area and higher 
time resolving power be conceived and a systematic attention 
be given to closer-by GRB sources. 

Despite this GRB is often consider ed as "unusual" 
JWatson et alJl2004l ISoderberg etai1 l2004). in our treatment 
we are able to explain its low gamma-ray luminosity in a 
natural way, giving a complete interpretation of all its spec- 
tral features. In agreement to what ha s been concluded by 
ISazonov. Lutovinov" & Sunvaev (2004a), it appears to us as 
a under-energetic GRB {Ed n , ~ 10 50 erg), well within the 
range of applicability of our theory, between 10 48 erg for GRB 



980425 jRuffini et al.ll2004bT) and 10 54 erg for GRB 991216 
jRuffini et alJ 

The precise knowledge we have here acquired on GRB 
031203 will help in clarifying the overall astrophysical sys- 
tem GRB 031203 - S N 20031w - the 2- 10 keV XMM and 
Chandra data (see e.g. Ruffi ni et all20 05c). 
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